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Abstract 

We investigate the dynamics of a relativistic electron in a strongly nonlinear plasma wave in 
terms of classical mechanics by taking into account the action of the radiative reaction force. The 
two limiting cases are considered. In the first case where the energy of the accelerated electrons 
is low, the electron makes many betatron oscillations during the acceleration. In the second case 
where the energy of the accelerated electrons is high, the betatron oscillation period is longer than 
the electron residence time in the accelerating phase. We show that the force of radiative friction 
can severely limit the rate of electron acceleration in a plasma accelerator. 
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1. Introduction 



New acceleration methods that provide a high rate of energy gain by particles are cur- 
rently being developed to generate charged particles with very high energies. One of the 
promising methods is the acceleration of charged particles in a plasma wave excited by a 
short intense laser pulse or a dense bunch of relativistic electrons . The longitudinal 
electric field in a plasma wave can reach huge strengths. Such a field is capable of accelerating 
electrons to very high energies. 

Recently, substantial progress has been made in generating short quasi-monoenergetic 
beams of ultrarelativistic electrons in laser plasma [5| . One of the models that describe 
the generation of a quasimonoenergetic beam of ultrarelativistic particles suggests that this 
generation is associated with the transition to a strongly nonlinear regime of laser pulse- 
plasma interaction. In this regime, the periodic plasma wave in the wake of a laser pulse 
transforms to a plasma cavity without electrons. Plasma electrons, along with external beam 
electrons, can be trapped by the cavity and accelerated to very high energies. 

It should be noted that, apart from a strong accelerating electric field, significant focus- 
ing transverse fields act on the electrons. These forces may be comparable in intensity to 
the accelerating force. The action of the focusing forces leads to the excitation of betatron 
electron oscillations across the acceleration direction. As a result of these oscillations, the 
relativistic electrons intensely emit electromagnetic waves jsl . Just as in conventional accel- 
erators, the losses through radiation can greatly reduce the electron acceleration efficiency 
in a plasma cavity. 

Since the properties of the electromagnetic radiation from an electron are determined 
mainly by the action of the focusing forces, when these properties are studied, the action of 
the accelerating longitudinal force may be disregarded, i.e., the electron dynamics may be 
considered in an infinitely long ion channel. In this case, the shape of the electron trajectory 
in the channel is nearly sinusoidal [9^ and is defined by the betatron frequency. 

^=^, (1) 

where Up = (47re^no/m)^^^ is the electron plasma frequency, 7 is the relativistic gamma 
factor of the electron, tiq is the unperturbed plasma density, e is the electron charge, and m 
is the electron mass. 
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The trajectory of an electron in an ion channel is similar to the helical trajectory of an 
electron in a constant magnetic field. The electron gyration in a magnetic field is defined 
by the cyclotron frequency u^. = eH/{mc) [11], where H is the magnetic field strength and 
c is the speed of light. As a result, the electromagnetic radiation spectrum of a relativistic 

in an ion channel is similar to the synchrotron 



radiation spectrum [ll| of an electron in a magnetic field. In the case of a relativistic 
transverse electron momentum, the spectrum is quasicontinuous. The synchrotron radiation 
spectrum is defined by a critical frequency that is a function of the betatron frequency for 
the radiation from the ion channel 
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fujJcr = -'J^hro^l'^/c ~ 5 X 10"^^7^no [cm~^] tq [/xm] keV, (2) 

where tq is the amplitude of the betatron electron oscillations in the channel. For frequencies 
much lower than the critical frequency, the radiated energy increases as ~ cj^/^, reaching its 
maximum near ~ 0.29uJcr, and decreases exponentially at > ujcr- 

By analogy with the synchrotron radiation of an electron in a magnetic field, which was 
studied in detail in jl2l], the effects related to the action of the radiative reaction force on 
an electron moving in an ion channel can be estimated. When the radiated photon energy 
becomes comparable to the electron energy ( mc^7 < Hucr), the electron dynamics becomes 
quantum one. The limiting electron energy at which the quantum effects become significant 
can be estimated using Eq. ([2]) 

^ ^ 4 mc' 
3 Hcroujp 

Strictly speaking, the quantum nature of the radiation can lead to a broadening of the 



electron orbit even at energies below the limiting value of (I3l) p^. However, these effects, 
which are significant for the long-term electron dynamics characteristic of conventional ac- 
celerators, are insignificant for the acceleration in plasma, because the interaction time is 
short. 

The space-time distribution of the electromagnetic fields in a plasma cavity are discussed 
in Section. 2. We derive the equations that describe the dynamics of a relativistic electron in 
a plasma cavity by taking into account the radiative reaction force. We obtain an estimate 
for the limiting energy of the accelerated electron that is limited by radiative losses. In 
Section 3, the electron dynamics in a plasma cavity is considered by taking into account the 
radiative reaction force in the case where the electron makes many betatron oscillations while 



flying through the cavity. We investigate the evolution of the electron energy as a result 
of radiative losses in an infinitely long ion channel. In Section 4, the electron dynamics is 
analyzed in the opposite limit of high energies, where the electron without the action of 
the radiative force makes no more than one betatron oscillation while flying through the 
cavity. In Section 5, we present the results of our numerical simulations of the acceleration 
of a relativistic electron bunch in a plasma cavity by taking into account the radiative 
reaction force using a hybrid particle- in-cell code (software package). The results obtained 
are discussed in Section 6. 



2. ELECTRON DYNAMICS WITH THE RADIATIVE REACTION FORCE 



As follows from our numerical simulations, the shape of the plasma cavity generated in 
the wake of a laser pulse is nearly spherical. The velocity of the plasma cavity, vi, is equal 
to the group velocity of the laser pulse and is close to the speed of light. The ions may be 
considered to be stationary, since the characteristic ion response time Tj ~ (47re^no/mj) 
for typical interaction parameters is much longer than -R/c, where rrij is the ion mass, R is 
the cavity radius, c is the speed of light. Thus, the electron density in the cavity is zero, 
while the ion density is hq. In this case, the potential distribution inside the cavity is [3] 

^ = -A^ = , (4) 

where ^ = x — vit, ip is the scalar potential, A is the vector potential, and the gauge 
if = —A^is written. In what follows, we use dimensionless quantities where the time, the 
velocity, the length, the electromagnetic field strength, and the electron density, n, are 
normalized to 1/ujp, c/up, majjp/\e\, and no , respectively. The laser pulse and the plasma 
cavity are assumed to propagate along the x axis. The space-time distribution of the quasi- 
static plasma fields inside the cavity is described by a linear function of the coordinates and 



time 



7| 



= e/2, Ey = -B, = y/A, E, = By = z/A. (5) 
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Let us first consider the electron acceleration in the cavity without the radiative reaction 
force. For the sake of simplicity, we assume that the electron trajectory lies in the z = 



plane and that the electron is accelerated along the x axis, making betatron oscillations 
along the y axis, with ^ Vy ^ ^i where px and Py are the longitudinal and transverse 
electron momenta, respectively. In this case, it is easy to find the trajectory of the electron 
acelerated in the cavity 

^ ^ + 2^' y^r,(^-^y co^ nd?j , (6) 

Here, the energy of the electron is assumed to be 70 at the time it enters the plasma cavity 
^ = —R, Tq is the initial electron deviation from the cavity axis, and 7; = (1 — vf)^^^"^ is 
the relativistic gamma factor of the laser pulse. Because of the electron acceleration, the 
betatron frequency Q = l/y^2rf is a slowly varying function of the time. The dependence of 
the electron energy on ^ is defined by 

7^7o + |(i?^-a. (7) 

The electron energy increment, 7 70 + 7f-R^/2, , is at a maximum at the cavity center 
= 0). It follows from Eq. ([6]) that the amplitude of the betatron electron oscillations 

decreases as the cavity center is approached. 

Making betatron oscillations, an ultrarelativistic electron emits electromagnetic waves 

and undergoes recoil as a result of the photon emission. The relativistic equations of motion 

for an electron in an electromagnetic field with the radiative reaction force are 11| 

^ = r'u, + ^,g\ (8) 
9' = -Q^^ku' - r'Fkiu^ + {Fkiv}) {F^^Um) u\ (9) 

where Fn^ Fik is the electromagnetic field tensor, is the 4-velocity of the electron, /i = 
2reCi;p/(3c) , r^, = 3-10^^ is the classical electron radius. The first term in Eq. ([8]), corresponds 
to the Lorentz force and the second term describes the action of the radiative reaction force. 

Under our assumptions, the focusing forces (the transverse component of the Lorentz 
force) make a major contribution to the energy losses through radiation, while the accelera- 
tion is attributable to the action of the longitudinal component of the Lorentz force. In this 
case, the approximate equations that describe the electron dynamics in the cavity are 

dp y 2 
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dy 

dt 7' 



(11) 




(12) 



The first two equations describe the betatron electron oscillations. When the force of radia- 
tive friction is disregarded (/i = 0) and for a slowly varying betatron frequency, the solution 
of the first two equations is ([6]). The first term on the right-hand side of the last equa- 
tion describes the action of the longitudinal component of the Lorentz force that provides 
electron acceleration, while the second term describes the radiative losses. 

Clearly, acceleration is possible as long as the accelerating force in the cavity is larger than 
the force of radiative friction. The limiting energy of the electron at which its acceleration 
is possible can be easily estimated from Eq. (fT2il 



Inequality (fT3|) is written in dimensional units to take into account the influence of the 
plasma density in explicit form. It follows from the derived inequality that the square of the 
limiting energy of the electron at which its acceleration is still possible is proportional to the 
cavity size and inversely proportional to the plasma density and the square of the distance 
to the cavity axis. The stronger the focusing of the beam, the higher the energies to which 
it can be accelerated in the cavity. 

Below, we will consider the two limiting cases. In the first case, the energy of the electron 
is relatively low and it makes many betatron oscillations while flying through the cavity. In 
the second case, at fairly high electron energy, the time of flight of the electron through the 
cavity is shorter than the betatron oscillation period. 

3. ELECTRON DYNAMICS AT A LARGE NUMBER OF BETATRON OSCIL- 
LATIONS 

Let us first consider the case where the electron makes many betatron oscillations while 
flying through one half of the cavity with an accelerating electric field. This condition limits 
the electron energy above. 




(13) 




(14) 
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FIG. 1: Transverse coordinate of the electron y (a) and its energy 7, (b) versus time for the motion 
in an ion channel with the parameters 70 = 5 • 10^, tq = 0.2 uq = 10^^^'^. Line 1 (solid) and line 2 
(dashed) correspond to the approximate solution ([25]) and the numerical solution of the equation 
of motion ([8|), respectively. 

For typical parameters of the interaction between an intense laser pulse and plasma 6| 
i? ~ 5, 7; ~ 10 condition ( fT4|) is satisfied for electrons with energies below 10 GeV. 

When the number of betatron oscillations is large, it is convenient to use the averaging 
method llSl. To this end, let us introduce a new variable. 



b exp ( i I Vldt ) ~ ^ ~ iVlp. 



(15) 



Substituting Eq. ( |T5l) in Eqs. (|TT i) -( fT2l) and averaging over the phase yields the following 
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reduced equations: 



(16, 



dt 2 V 272^1 2 



R-7^2]-^\bf7'. (17) 



It should be noted that although rapidly oscillating component of 7 is small the derivative of 
it with respect to time can be large. So to derive Eq. (fT6ll we use Eq. (fT2l) before averaging 
in order to exclude d'y/dt. 

Assuming that s = \bf, we can rewrite Eqs. ( IT6l) and (fT7|) as 

'^■5 ^^ 2 /1o^ 

? = -f^-AVT^7^ (19) 

where ro = a/2 (y^) = 2a/s is the amplitude of the betatron electron oscillations in the 
cavity. In the absence of an accelerating force, the equations take the form 

These equations describe the electron dynamics in an infinitely long ion channel where the 
longitudinal component of the Lorentz force may be disregarded. The system of equations 
( IT6l) and (|T71) has the integral of motion 

= const = (22) 

where the integration constant can be found from the initial conditions. It follows from 
Eq. ( 122|) that the betatron oscillation amplitude slowly decreasing as ?/ oc 7^/^ in this ap- 
proximation. Thus, the radiated energy is pumped out from the longitudinal component as 
well as from transverse component of the electron energy. Solving Eqs. (fT6!) and (fT7|l yields 

^ = (l + 5/iro?7o/32)4/5- ^^^^ 
As an example, let us consider the radiative deceleration of an electron with the initial 
parameters 70 = 5 ■ 10"^, ro = 0.2 in an ion channel with the density uq = 10^^~^. As follows 
from Fig. 1, the constructed solution is in good agreement with the results of our numerical 
integration of the equation of motion ([8]) . 
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Let us now consider the electron dynamics in a plasma cavity by taking into account 
the accelerating force and the radiative reaction force. Let us assume,that the betatron 
oscillation amplitude is approximately constant with time. This assumption is valid when 
the energy increment of the electron as a result of its acceleration is considerably smaller 
than its initial energy (A7/7 ^ ifR^ / {^l) ^ !)• Integrating Eq. (IT^ with the condition 
{y'^) = rl/2 = const yields 

^ 4a At' jaQ + SBt' {aQ , . 

^ fir^o-ff A% «) + «) ' 

where ^ = —R + t/{2'-)f), Ai{x) and Bi{x) are the Airy functions iJ] , Ai'{x) and Bi'{x) 
are the derivatives of the Airy functions, and a = —'yf^^rQ^^jj^^'^2~'^^^, 6 is the integration 
constant, which can be found from the initial condition 7 = —R) = 70. It should be noted 
that the factor a does not depend on the electron energy. Moreover, for typical laser-plasma 
interaction parameters, a.^ <^ 1, and the asymptotics of the Airy function may be used for 
low values of the argument. Expanding separately the numerator and the denominator in 
Eq. (1241) in terms of the small parameter a^, we obtain the following expression for the 
electron energy: 

^ 70 (32 + fiRt'r',) + A^ft {4jfR - t) 

32 + /itr2 (470 + ti?) ■ ^ ' 

Let us consider the electron dynamics in a plasma cavity with the parameters 70 = 4-10^, 
7i = 10, i? = 15, ro = -R/4, riQ = 10^^~^. Figure 2 presents the time dependences of the 
electron density obtained by integrating the equation of motion (jHj) with and without the 
radiative reaction force. At t = 3000, the electron reaches the cavity center, where the 
accelerating force vanishes. As we see from Fig. 2, when the radiative losses are disregarded, 
the electron energy reaches its maximum at the cavity center. Including the radiative losses 
leads to a significant reduction in the maximum electron energy. It follows from Fig. 2 
that the approximate solution (125!) is in good agreement with the results of our numerical 
integration of the equation of motion. 



4. ELECTRON DYNAMICS IN THE ABSENCE OF BETATRON OSCILLA- 
TIONS 

In the opposite limit of high electron energies (7 > R'^'yf /n"^) without the action of the 
radiative force, the electron makes no more than one betatron oscillation while flying through 
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FIG. 2: Electron energy versus time for the motion in a plasma channel with the parameters 
70 = 4 • 10"^, 7i = 10, R = 15c/ujp, tq = i?/4, no = 10^^~^. The solid, dashed, and dotted lines 
correspond to the approximate solution (j25p . the numerical solution of the equation of motion ([8]), 
and the numerical solution of the equation of motion ([8]) without the action of radiative friction 
{fj, = 0), respectively. 

the accelerating half of the cavity. Let us assume that the accelerating force in the cavity 
is larger than or comparable to the radiative reaction force in the longitudinal direction, 
7^ < 2R/ {^Tq) . Since the energy increment during the acceleration at high electron energies 
is relatively small, the electron energy in the equations for the transverse coordinate ( JTTl) and 
f[T^ may be considered as a constant. Thus, for the transverse coordinate of the electron, 
we derive the expression 

y it) ^ ro cos Vlt. (26) 
Substituting this expression in Eq. (fT2!l yields an equation for the electron energy: 

'^ = \(r-1\- ^^2 2 C0S2 nt. (27) 

dt 2 V 272^ 4 ' ° ^ ^ 

For the case where fit ^ 1 , the solution of this equation was found in the previous section 
and is represented by Eqs. ( !24l) and (!25|) . in which the substitution Tq — > 2r1 should be 
made. 

As the initial electron energy increases, the accelerating force can become negligible 
compared to the force of radiative friction (7^ ^ 2Rjj,~^rQ^). In this case, the energy losses 
through radiation can be significant. The betatron oscillation amplitude increases as the 
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electron energy decreases as a result of the generation of electromagnetic radiation and the 
time dependence of the transverse coordinate should be taken into account. By disregarding 
the accelerating force, we can reduce Eqs. (ITT1) - (IT^ to the form 

% - (''^ 

It should be noted that in the absence of an accelerating force, the terms in Eq. (125]) that 
describe the action of the nonlinear force of friction cancel each other out. 

In the zeroth approximation, we will assume that the transverse coordinate does not 
depend on the time. The time evolution of the electron energy in the cavity is then defined 
by Eq. ( l23l) . In the next approximation, we will take into account the time dependence 
of the electron energy in Eq. (1281) . As a result, we obtain the following equation for the 
transverse coordinate: 



Its solution can be expressed in terms of the Airy functions: 

y = Trro [At (r) Bt' (e) - Bi (r) Ai' {e)] , (31) 

r = £(l + /5t), (32) 

where e = —2jj,^'^/'^'yQ^rQ^^^, (3 = n^y^r'^/A. Substituting the solution obtained in Eq. (12^ 
yields an equation for the electron energy. The solution of this equation can also be expressed 
in terms of the Airy functions: 



-- — = — {t\Ai (r) Bi' {e) - Bi (r) Ai' [e)? + [At! (r) Bi' {e) - Bi' (r) Ai' {e)f - ^\ . 
7 7o 1q£ ^ TT^J 

(33) 

In the limit /5t <^ 1, the solution obtained transforms to Eq. (123|) . 

As an example, let us consider the motion of an electron with the initial energy 70 = 4-10^ 
and the initial deviation from the axis tq = 0.1. The ambient plasma density is = 10^^~^. 
Figure 3 shows the dynamics of the transverse coordinate of the electron and its energy 
for the motion in a cavity with the reaction force. As follows from Fig. 3, the derived 
expressions (I3T!) and (133|) are in good agreement with the results of our numerical integration 
of the equations of motion. 

11 




1000 



1000 



2000 3000 



(Opt 



b) 




2000 3000 



(Opt 



FIG. 3: Transverse coordinate of the electron (a) and its energy (b) versus time for the motion in 
a plasma cavity with the parameters 70 = 5 • 10^, tq = 0.2 rig = 10^^~^. The solid, dashed, and 
dotted lines correspond to the approximate solution (j31|) and (j33|) . the numerical solution of the 
equation of motion (l8|), and the approximate solution ()23p and y =const, respectively. 



5. RESULTS OF NUMERICAL SIMULATIONS 



The action of radiative friction on the external beam of relativistic electrons accelerated 
in a plasma cavity was numerically simulated using a two-dimensional relativistic hybrid 
particle- in-cell code for a cylindrical geometry 18|. The quasi-static approximation (the 



plasma response is assumed to vary slowly during the characteristic evolution time of the 
laser pulse) is used to speed up the computation. The code includes the action of the force 
of radiative friction on the electron dynamics. 
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In our numerical simulations, we assumed that the laser pulse was circularly polarized 
and had a Gaussian profile, a = aoexp {—r'^/rf — C,'^/Lf), and the wavelength A = 0.82 m. 
The laser pulse parameters are r; = 5, Li = 2, oq = 10. The pulse propagates in a plasma 
with the density uq = 10^^^ cm~^. The beam parameters are the following: the relativistic 
gamma factor of the electrons is 7 = 10^, the beam radius is = 2, and the beam electron 
density is Ub = 10"^'' cm~'^. 

Figure 4a shows the distribution of the electron plasma and beam density. The beam 
electrons are accelerated in a plasma cavity with the radiative reaction force. The transverse 
forces in the cavity focus the beam. As follows from Fig. 4a, the beam electrons reach the 
cavity center approximately in half the betatron period, which matches the estimate of the 
time it takes for the electron to reach the cavity center. 

The beam electron distribution function at the time the leading edge of the beam reaches 
the cavity center is shown in Fig. 4b (line 1). At the beginning of the interaction, all 
of the beam electrons had the same energy (line 2 in Fig. 4b). As follows from Fig. 4b, 
the radiative reaction force decelerates the beam electrons with a large betatron oscillation 
amplitude, while the electrons with a small amplitude are accelerated by the longitudinal 
electric field. According to our numerical simulations, the energy of the electrons with a 
small betatron oscillation amplitude increases by 1 GeV, which is close to the estimate that 
follows from Eq. ([7]). The force of radiative friction for the electrons with the maximum 
oscillation amplitude (ro — 2) exceeds the accelerating force; as a result, these electrons 
are decelerated, losing 2 GeV of their energy (see Fig. 4b). A similar estimate (2.5 GeV) 
for the energy lost by the electron through radiation is obtained from Eq. (J25l) . where it is 
considered that t ~ -n^jT?). 

6. CONCLUSIONS 

Figure 5 shows the various regimes of electron dynamics in a plasma cavity with the 
radiative reaction force in the 70 (initial electron energy)- tq (initial deviation) plane. As 
follows from Fig. 5 and inequality f|T3|) . the force of radiative friction can restrict significantly 
the electron acceleration in plasma. The radiative friction can be reduced by reducing 
the beam radius and the plasma density. For the plasma density = 10^^"'^, electron 
acceleration is not possible at energies above 50 GeV and at a radius of the accelerated 
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FIG. 4: (a) Plasma electron and electron beam density. Darker colors correspond to higher 
electron densities. The location of the laser pulse is indicated by the dotted line. The coordinates 
are presented in units of c/oop. (b) The beam electron energy distribution function: at the initial 
time (dashed line 1) and at the time the leading edge of the beam reached the cavity center (solid 
line 2). 



bunch To > 10 m, which is typical of present-day accelerators [19|. However, the projected 
accelerators (SLAC, TESLA) are expected to provide beam focusing up to tens or hundreds 



of nanometers 
TeV. 
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2l| . In this case, electron acceleration is possible up to several tens of 



Since the radiated photon energy increases quadratically with electron energy, these ener- 
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FIG. 5: Various regimes of electron dynamics in a plasma cavity with the radiative reaction 
force in the 70 (initial electron energy)-ro (initial deviation) plane. Straight line 1 corresponds to 
the boundary at which the accelerating force is equal to the force of radiative friction. Straight 
line 2 corresponds to the boundary at which the time of flight through the cavity is equal to 
the betatron oscillation period. Straight line 3 corresponds to the boundary that separates the 
quantum domain of radiation from the classical one. The interaction parameters are 7^ = 10, 
R = Sc/wp, no = 10^^-^. 

gies can become equal for ultrarelativistic electrons. In this case, the radiation and the elec- 
tron dynamics itself are of a quantum nature. As follows from Fig. 4, for uq = lO^^cm"'^ and 
tq > 10, the quantum effects become significant for electrons with energies above 500 GeV. 
Thus, the quantum theory should be used at high electron energies. In recent years, im- 
pressing results on the acceleration of electrons in plasma by relativistic electron bunches 

le ion-focusing regime where 



23J. In this case, the plasma 



have been obtained |22] . The acceleration is most efficient in t^ 
the bunch density is higher than the ambient plasma density 
electrons are expelled almost completely from the region around the bunch by its large neg- 
ative charge. The electromagnetic field distribution in this region is identical to the field 
distribution inside the plasma cavity generated by a short laser pulse p/7|]. Thus, the re- 
sults obtained are also valid for describing the dynamics of the electrons accelerated by a 
relativistic electron beam in plasma. 
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